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We treat a fault as an array of asperities with a prescribed statistical 
distribution of strengths. When an asperity fails the stress on the failed 
asperity is transferred to one or more adjacent asperities. For a linear 
array the stress is transferred to a single adjacent asperity and for a 
two-dimensional array to three adjacent asperities. Using a renormalization 
group (RG) method to extrapolate to an arbitrarily Large scale we show that 
the solutions bifurcate at a critical applied stress. At stresses less than 
the criticaL stress virtually no asperities faiL on a large scale and the 
fault is locked. At the critical stress the soLution bifurcates and asperity 
failure cascades away from the nucleus of failure; we interpret this 
catastrophic failure as an earthquake and it corresponds to the transition 
from stick to slip behavior on the fault. Thus the stick-slip behavior of 
most faults can be attributed to the distribution of asperities on the 
fault. Our results explain why stick-slip behavior on faults is commonly 
observed rather than stable sliding, they explain why the observed level of 
seismicity on a locked fault is very small; and they explain why the stress 
on a fault is less than that predicted by a standard value of the coefficient 


of friction. 


Introduction 


This paper presents a simple friction model to provide an understanding 
of the stick-slip behavior of faults. A large fraction of active faults 
appear to behave in a stick-slip rather than a stable sliding manner. For a 
fault system such as the San Andreas, which has well over 1000 km of surface 
exposure, approximately 55 mm/yr of relative motion, and has been extensively 
studied (Turcotte, 1977), it is possible to differentiate between these two 
alternate behaviors. The San Andreas can H*: subdivided into three principal 
sections based on the type of behavior observed: northern and southern 

locked sections and a central creeping section. There is strong evidence 
that the northern and southern locked sections, the sites of the 1906 and 
1857 earthquakes, behave in a stick-slip fashion. There has been no 
significant seismicity on these sections of the San Andreas in the period 
since adequate instrumentation has been available. In the cent. one, 

aseismic creep has been observed along with many small and moderate size 
earthquakes, suggesting a stable sliding type of behavior. 

Since repetitive sequences of earthquakes are observed to occur on 
active fault systems it is appropriate to treat faults as approximately 
planar surfaces with a coefficient of friction. Using frequency magnitude 
and moment magnitude relationships Aki (1981) has shown that the fractal 
dimension of a fault is D = 3b/c, where b is the slope of the log frequency 
magnitude relation and c is the slope of the log moment magnitude relation. 
For c = 1.5 (Hanks and Kanamori, 1979) and b = 1, the fractal dimension is 2, 
the same as the topological dimension of a plane (Aki, 1981). For regions 
where the b value is near 1 the planar approximation should therefore be 
quite good. Since most earthquakes occur on preexisting faults it is not 
appropriate to model earthquakes on active faults as the fracture of pristine 


rock. While surficial fault traces have been observed to go through 
previously unfractured rock, there is much evidence that earthquakes 
repeatedly occur on the same fault surface or zone of surfaces (Sieh,, 1 9 7 & ) . 

We will therefore model a fault system which can contain fault trace 
offsets, multiple fault traces and varying amounts of curvature or roughness 
as an array of asperities with a statistical distribution of strengths. 

There is evidence that large asperities, or barriers, can control the 
propagation of an earthquake along a fault. 

Several authors have previously modelled the stick-slip behavior of 
faults in terms of frictional effects. Weertmann (1979) modelled the 
instability in terms of a frictional stress on a fault that decreases with 
increasing slip velocity. Stuart and Mavko (1979) modelled the instability 
in terms of a strain-softening constitutive relation for the fault zone. 

In this paper we present an alternative hypothesis to explain the stick-slip 
behavior of faults. 

In terms of mathematical modelling an earthquake is clearly a 
catastrophic change in the behavior of the system. Recently, renormalization 
group (RG) techniques have been successfully applied to models that exhibit 
catastrophic behavior of the type found in natural systems (Wilson and Kogut , 
1974; Fisher, 1974). A classical example is that of a system undergoing a 
phase transition. 

Renormalization group (RG) techniques have been used by Madden (1981) to 
relate the macroscopic electrical conductivity and fracture of rocks to the 
microcrack population, and by Allegre et si, (1982) to study the coalescence 
of fractures. Newman and Knopoff (1982, 1983) have also studied the 
coalescence of fractures and while they use the term renormalization in their 
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work no rescaling is done so their approach is substantially different from 
the usual RG methods. 


Formulation of the Problem 

In this paper we model a fault as an array of asperities with a 
statistical distribution of strengths. We will consider both linear and 
two-dimensional arrays of asperities. In order to illustrate the approach we 
will first consider a linear array as illustrated in Figure 1. This model 
should be appropriate for large scale asperities (barriers) on a long fault. 
The fault is broken into n elements of length <5x and each element is assigned 
an asperity failure strength f; the asperity will fail when the stress on 
the asperity reaches this value. The asperities have a distribution of 
strengths which will be specified by a statistical distribution function. 

When a stress is applied to the fault all asperities with a failure 
strength f < will fail. 

We will divide the linear array of n asperities into n/2 cells, each 
containing two asperities as illustrated in Figure 1. When one asperity in a 
cell fails we assume that the stress on that asperity is transferred to the 
other asperity in the cell. This is an essential feature of our model and is 
equivalent to the transfer of stress to adjacent regions when a crack is 
introduced into an elastic solid. The cell size is a measure of the distance 
over which the stress is redistributed after an asperity failure. The reason 
we assume that there are two asperities in a cell is that the stress in the 
failed region is applied over a length which is of the order of the length of 
the failed region. 

A cell may contain two broken, two unbroken, or a broken and an unbroken 
asperity. If a cell contains a broken and an unbroken asperity the strength 
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a ^ of the unbroken asperity must be greater than 2 a if it is to survive 
the transfer of stress a from the broken asperity. We assume that a ceil 
fails only when both asperities in the cell fail. The basis of the RG 
approach is that after the first renormalization, each first order cell is 
now treated as a second order (r =* 2) asperity and pairs of second order 
asperities form second order cells. This is illustrated in Figure 1. The 
statistics of failure of the second order asperities and cells is the same as 
the first order asperities and cells. The process is repeated by iteration 
to infinite order. It should be emphasized that the same statistical 
distribution of asperity strenghts is applied to the higher order cells even 
though asperities are destroyed at each of the lower levels of calculation. 
This is an essential feature of the RG method. 

The process of stress transfer and induced failure tends to increase the 
lengths of segments of broken asperities. As the applied stress is 
increased, a value is reached at which failure of an infinite length of 
asperities will occur and the behavior changes catastrophically from stick to 
slip. The stress at which this change occurs is equivalent to the 
temperature at which a phase transition occurs. The statistical distribution 
of energies in a solid, liquid, or gas is equivalent to the statistical 
distribution of asperity strengths in our model. The utilization of the RG 
method allows us to study the development of failed segments as the 
characteristic lengths of the failed segments increase with increased applied 
stress , 

We first consider the distribution of asperity strengths. Clearly a 
wide variety of strengths on a wide variety of scales must exist on any real 
fault. For example, fault bends and offsets correspond to strong 
asperities. However, data on actual distributions of asiperity strengths are 
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not available. It is possible that studies of the type given in this paper 
may allow asperity distributions to be inferred from such seiflmic 
observatiors as the dependence of earthquake frequency on magnitude. In the 
absence of applicable data we assume a quadratic Weibull distribution for the 
probability P a that the failure strength o f of an asperity is less than 
the stress ao 


P * 1 
a 


-(ax) 

e 


2 


(i) 


where 

x = a/a 

o 

and a 0 is a reference asperity strength. Wiebull distributions are often 
used to represent a statistical distribution of failure strengths (Harlow and 
Phoenix, 1982). It should be emphasized that our approach can be applied to 
any continuous distribution of asperity strengths. The probability 

2 

P x » 1 - e" X (2) 

that a f a is shown in Figure 2a as a function of a/o 0 . Ten percent of 
the asperities have failed when o/a Q - 0.32, fifty percent of the 
asperities have failed when a/a 0 = 0.83, and ninety per cent of the 

asperities have failed when a/o 0 = 1.52. Since the relation between Pi 

and a is invertible, can be used as a measure of the applied stress. 

The probability that failure will occur at the applied stress a/a 0 is given 
by dPi/dx and is shown in Figure 2b. The probability that Of = a is zero 
at zero stress and increases to a maximum at a = 0.71 a 0 . The mean 
strength of an asperity is a = (/tt/ 2 )a Q = 0.8862 a Q . 

An essential feature of our model is the transfer of stress from a 
failed asperity to its nearest neighbors. Without this transfer of stress 
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the behavior of the system is simple and uninteresting. Since strong 
asperities will not break until large stresses are applied they can block the 
propagation of broken segments and there is no change from stick to slip 
behavior on the fault. It is the transfer of stress from broken asperities 
onto the remaining unbroken asperities that leads to catastrophic behavior at 
an applied stress that is less then the average strength of the asperities. 

It is clear that stress transfer to the adjacent unbroken sections of a fault 
will occur on a real fault. 

In order to quantify the failure of asperities due to the transfer of 
stress we introduce the conditional probability P a> b that failure will 
occur when a stress (a-b)a is transferred to an unbroken asperity supporting 
a stress bo, so that the final stress on the asperity is ao . This 

i 

conditional probability is given by 


P 


a,b 



(3) 


with . 

2 

P * 1 - (1 - Pi) 3 (4) 

a 1 

for the probability function given in (1). 

In principle this problem could be solved without the use of the RG 
technique. However, the range of scales that could be studied is quite 
limited even with the largest computers available. Although we will utilize 
the renormalization group method in the standard manner, it should be 
recognized that the approach is semi-empirical and has been principally 
justified by its success in solving a variety of fundamental unsolved 
problems in physics. These problems fall in a broad class in which a 
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continuous system on the microscopic scale exhibits catastrophic behavior on 
the macroscopic scale. We argue that the stick-slip behavior of faults falls 
in this general classification of physical problems. 

We will first illustrate the application of the renormalization group 
technique to the failure of a linear array of asperities using a basic cell 
composed of two •'asperities” and the probability distribution given in (1). 
Note that the (r+l)th order "asperities” which result from r iterations of 
the RG transformation contain 2 r actual first order asperities. The first 
three renormalizations are illustrated in Figure 1. For a cell containing 
two asperities which are either broken or unbroken, four states are possible: 
1) [bb], 2) [bu], 3) [ub], and 4) [uu] , where b represents a broken asperity 

and u represents an unbroken asperity. Note that states 2 and 3 are 

* 

equivalent and can be combined into a single state with a multiplicity of 2. 
The probabilities for each of these states neglecting any interactions 
between asperities is given in Table la. 

Next it is necessary to consider the influence of a broken asperity on 
an adjacent unbroken asperity. We use the conditional probability ?2 ,1 
that an unbroken asperity already supporting a stress a will fail when an 
additional stress o is transferred to it from an adjacent broken asperity. 
Including the effects of such induced failures leads to the probabilities 
given in Table lb for each of the cell states. Since "the conditional 
probability from (3) is given by 


P 


2,1 



(5) 


the probabilities with stress interactions in Table lb can be expressed in 
terms of Pi and • 
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We must prescribe a condition for determining whether a rth order cell 
is broken or unbroken. We assume that a rth order cell is broken only if 
both "asperities" in the cell are broken. Under this condition the 
probability that a first order cell is brok^M Pj(2) 

is given by 

P x (2 * - Z? X V 2 - P 1 2 , (6) 

and substitution of P 2 from (4) gives 

p/ 2) * 2? x [l - (1 - P x ) 4 ] - V* . (7) 

For higher order cells (7) is used as an iteration equation to determine 
Pj(r+1) from where r is the order of the cell being 

considered. Implicit in the RG method is the assumption that the probability 
distribution applicable to a first order, r+1, cell is also applicable to 

i 

higher order cells, This assumption is an essential feature of the method 
and is clearly an approximation. The general form of (7) is 

Pl (r+1) = 2 P 1 <r) tl - (1 - - (P, (r) ) 2 • (8) 

The dependence of p^Cr+l) on Pi (r > is given in Figure 3. The points 
0 and 1 are stable fixed points of the system. The straight line 
corresponding to Pj(t+1) = Pj(r) is also included in Figure 3. The 
iterative relation crosses this straight line at Pj( r ) = P* s 0.2063. We 
will show that P* is an unstable fixed point that separates the region of 
stick behavior from the region of slip behavior. 

The RG iteration can be performed graphically using Figure 3. For 
example, we take Pi = 0.6 and from (8) find Pi^2) = 0.8093. This cell 
behavior at order 1 now becomes the asperity behavior at order 2. To do this 
graphically a horizontal line is extended to the line pj(r+l) = P^ r ) 


~T 


M 

to reflect the total cell behavior at order 1 into the asperity behavior at 
order 2. Thus the probability of cell failure at order 2 i3 Pj < 3) - 
0.9615. This procedure is repeated to give P[(4) * 0.9985, etc., and the 
probability of failure rapidly approaches unity as the order is increased. 

On the other hand, if we take Pj ** 0.1 we find Pj(2) * 0.05-K78, Pi^) 

* 0.02184, Pj/4) ■ 0.00322, etc., and the probability of failure 
decreases towards zero as the order is increased. If P^ > P* failure occurs 
for infinite length scales and slip behavior results. If P^ < P* the 
behavior is stable and failure occurs only on the smallest scales. 

Bifurcation of the solution occurs at Pj * P* * 0.2063 and the critical 
stress leading to failure is a* « 0,4.807 a 0 from (l). The dependence of 

P X < r) on r for several values of P^ is given in Figure 4. The 

bifurcation of the soLution at P^ = P* » 0.2063 is clearLy illustrated. Note 

that the value of the critical stress is considerably less than the value of 

the mean strength of an asperity a = 0.8862 a 0 . 

The stable behavior of the system at a < a* can be characterized by a 
correlation length L which measures the maximum length over which failure 
occurs for Pi < P*. The rapid increase of L as the threshold is approached 
from below is described by a power low 

LMP*- P x r V , (9a) 

or equivalently 

L a- (a* - a) V , (9b) 

where v is the correlation length exponent (Wilson and Kogut, 1974). 

According to this result the magnitude of precursory seismicity would be 
expected to increase as the critical stress on the fault is approached. The 
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onset of catastrophic behavior at o * 0 * corresponds to the divergence of the 
correlation length I#. 

The correlation length exponent v is easily obtained from the RG 
transformation (8). Given the dependence of Pj_(r+1) on pj(r) f the 
slope of the curve in Figure 3 at Pj ■ P* is given by 
(r+1) 

* ( 10 ) 


dP 


A 5 


1 


dP 


(r) 


p » P* 
1 


As long as (P* - P^ r ^) « 1, a linear approximation to (8) is valid, and 
P* - P, (r * 1 ^ 

A * L__ . (11) 


P* - P 


(r) 


1 


It then follows that (Wilson and Kogut, 1974) 

A * b 1/v , (12) 

where b is the linear rescaling factor. For the b = 2 RG transformation that 
led to (8) we obtain A * 1.6189, so that v a 1.4388. 

So far we have considered only a linear array of asperities. We will 
next consider a two-dimensional array of asperities distributed uniformly on 
a planar fault as illustrated in Figure 5. We will divide the 
two-dimensional array of n asperities into n/4 cells each containing four 
asperities. The failure of individual asperities will be treated in the same 
way as in the linear case and (1) is assumed to be applicable. When one or 
more asperities in a cell fail we assume that the stress on those asperities 
is transferred equally to the remaining asperities in the cell. That is, if 
one asperity fails the stress on the three remaining asperities is 4a/3. We 
choose four asperities in a cell so that the stress in the failed region is 
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applied over a Length which is of the order of the Length of the failed 
region. We again assume that the ceil fails when all asperities in the cell 
fail ♦ 


A second order ceil is composed of four first order cells or second 
order asperities and, therefore, sixteen primary asperities as illustrated in 
Figure 5. The statistics of failure of the second order asperities and ceLls 
is the same as the firsf ~,rder asperities and cells. Again, the process is 

repeated by iteration to infinite order. The RG transformation thus 

constructed corresponds to a linear rescaling factor b 8 2 on a two- 

dimensional array. This case is considerably more complex than the linear 

array example considered above. Following the same procedure illustrated in 
Table l and using the definition of the conditional probability we find that 
the probaV ; i l ty that a cell fails is given by 


P l (2> " P l 4 + 4 P i 3(1 " P 1 )P 4,1 + ^i 2 * 1 " p i )2 t p 2 2 l 
+ 2 P 2,l (1 - P 2,1> P 4,2> + 4P 1 (1 - P 1 ,3 1 P 4/3 ! 1 


+ 3 P 


4/3, 1 (1 " P 4/3,l )P 4,4/3 * 3P 4/3,1 (1 “ P 4/3 ,1^ [P 2,4/3 


+ 2 P 2,4/3 (1 " P 2,4/3 /Jr 4,2 J 




(13) 


and introducing (3) we obtain 


P l <2> “ P l 4 + 4P 1 3(P 4 " P l } * 6P 1 2(P 2 ~ P l )2 + 12P 1 2(P 2 “ P 1 )(P 4 " P 2 ) 
+ 4P l (P 4/3 - P l )3 + 12P l (P 4/3 ' P 1 ,2<P 4 * P 4/3> 


+ l2P l (P 4/3 ' V (P 2 ' P 4/3 > 

+ 24P l< P 4/3 - P 1 )(P 2 " P 4/3 >(P 4 - V ‘ 

The dependence of p^(r + D on P]^( r ) shown in Figure 6 follows from 
introducing P a from (4) and using (14) as an iterative relation. 


( 14 ) 
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The general behavior of this two-dimensional case is the same as that of 
the linear example considered above. Again, an S-shaped curve is generated. 
The points 0 and 1 are stable fixed points. The crossing at P* “ 0.1707 
separates stick from slip behavior. From (1) the bifurcation of the solution 
occurs at a* ■ 0.43/7 o 0 , Thi'j is just about one-half the mean strength of 
the asperities a * 0.8862 o 0 <. We also find that A * 2.357 ; the correlation 
length exponent v * 0.8084 follows from (12) with b * 2. The quantitative 
differences between these results and those for the linear array, as listed 
in Table 2, illustrate the effect of the physical dimensionality on the 
critical behavior of the system. Note the decrease in both the critical 
probability P* and the correlation length exponent v with increasing D. 
Simpler percolation models exhibit the same trend when D is increased from 
two to three (Stauffer, 1979). 


Conclusions 

We have shown that a statistical distribution of asperity strengths 
leads to stick-slip behavior of a fault. The transfer of stress from failed 
sections of the fault to adjacent locked sections is an essential feature of 
our model. We have used the RG approach to obtain the behavior of the model 
as a function of applied stress, The main result is the existence of a 
finite critical stress below which the fault breaks are always bounded in 
their growth. 1 

We also find that the value of the critical stress is considerably 
smaller than the mean strength of the asperities. This may explain the low 
stress levels associated with displacements on the San Andreas fault. 
Laboratory studies of friction generally result in a coefficient of friction 
near 0.6 (Byerlee, 1978). However, the low measured heat flow adjacent to 



1 
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the San Andreas fauLt is strong evidence that the equivalent coefficient of 
friction on the fault is about a factor of two less than the laboratory value 
(Turcotte et al., 1980). We find that the critical stress on a fault is 
about a factor of two less than the mean strength of asperities. A similar 
result has been found for fibrous materials (Harlow and Phoenix, 1982). 

Clearly the model considered in this paper is based on a number of 
simplifications. These include: 

1) The assumption that the applied stress a on the fault is a constant 
is not a good approximation. The stress on an actual fault will have spacial 
variations. Stress concentrations are expected to occur at the edges of 
locked sections. These concentrations would be expected to initiate 
catastrophic failures on a fault. However, significant levels of stress will 
have accumulated on all sections of the fault and our model shows how a 
failure, once nucleated, can spread due to the transfer of stress. 

2) The form of the asperity strength distribution given in (1) is 
arbitrary. However, the RG method given here is applicable to any continuous 
distribution of asperity strengths. 

3) Failure on actual faults does not extend to infinity. There is ample 
observational evidence that strong asperities (barriers) can block the 
propagation of a zone of failure. This limiting behavior can be included in 
our analysis by utilizing a more complex relationship for the distribution of 
asperity strengths. 

4) The assignment of the same scale <5x to all asperities is a poor 
approximation. Clearly large asperities may have larger physical dimensions 
than smaller asperities. On the smallest scale asperities may have atomic 
dimensions whereas major barriers such as bends in a fault may have 
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dimensions of tens of kilometers. However, this additional complication 
would not be expected to affect the qualitative predicted behavior. 

We would argue that our results make physical sense. Once a broken 
patch on a fault starts to grow the transfer of stress is sufficient to break 
adjacent asperities as long as the distribution of strong asperities is 
sufficiently small. Obviously real fault breaks have finite lengths. We 
argue that a fault break Is terminated when a very strong asperity is 
reached. This is confirmed by the observation that long fault breaks are 
often terminated by bends or offsets in the fault trace. The scale upon 
which large asperities are distributed will determine the length of fault 
breaks and the behavior of the fault. 
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a) Probabilities without stress interactions 
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b) Probabilities with stress interactions 
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2P P - P 
1 1*2 *1 


1 + Pf - 2P l P 2 


c) Probability of failure applied to the next order cell 


Table 1. Renormalization group applied to a basic cell composed of two 
asperities. The probability of the failure of an individual asperity is Pj 
and the probability of an induced failure is P 2 . 
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D 

p* 

o*/o 0 

V 

1 

0.2063 

0.4807 

1.4388 

2 

0.1707 

0.4327 

0.8084 


Table 2. Renormalization group results for the critical probability P*, the 
critical stress a*/cf 0 , and the correlation length exponent v in one and two 
dimens ions . 


Figure Captions 

Figure l. Illustration of the statistical distribution of asperity 

strengths Of for a linear array; o 0 is a reference asperity 


21 


strength. An asperity is assigned to each unit length 5x. Also 
shown are the cell sizes for orders r = 1 to 4. 

Figure 2. (a) Dependence of the probability P^ that failure of an asperity 

will have occurred on the normalized stress o/o 0 . (b) 

Dependence of the probability that failure will occur at the 
normalized stress o/o Q . This is the change in the probability 
of failure SPj when there is a change in the normalized stress 
6 (a/a 0 ) . 


Figure 3. Dependence of the probability of failure for the r+1 cell 
Pf(r + i) on the probability of failure of the r cell Pj^ r ) 
for cells containing two asperities with a quadratic Weibull 
distribution of strengths. The procedure described in the text 
for determining the probability of cell failure for successive 
iterations is illustrated for P^ = 0.6, 0.1 . The critical 
probability of failure P* gives the bifurcation point for 
catastrophic failure of the system. If Pj^ < P* the solution 
iterates to Pj® = 0 and no failure occurs. If > P* the 
solution iterates to Pj® = 1 and the system” has failed. 

Figure 4. Dependence of the probability of failure P^ r ) on the order r 

for several values of Pj , The bifurcation of the solution at P^ 
P* = 0.2063 is clearly illustrated. 

Figure 5. Illustration of the two-dimensional array of asperities with 

four asperities per cell. Second (2), third (3), and fourth (4) 


order cells are also shown. 
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Figure 6. Dependence of the probability of failure for the r+1 cell 
P^(r-fl) on the probability of failure of the r cell Pj(r) 
for cells containing four asperities with a quadratic Weibull 
distribution of strengths. The critical probability of failure P* 
gives the bifurcation point for the system. If P^ < P* the 
solution iterates to Pi® = 0 and no failure occurs. If Pi > P* 
the solution iterates to Pi" = 1 and the system has failed. 
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